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Abstract 

We derive an exact equation for the critical velocity for He II superfluidity 
motion in channels of width < 10~ 5 m for which no quantitatively satisfactory 
prediction exists prior to this work. 

*This e-print version (2) contains fuller list of the earlier reports on the new superfluid theory. 



In channels of width d > 10 5 m, the critical velocity for He II superfluidity motion 
can be satisfactorily described by[TJ [2] 



y m a 

a prediction based on the QCE (quantum confinement effect) superfluidity mecha- 
nism. In ([T]), a (~ 3.6 x 10~ 10 m) is the average interatomic separation distance of 
He II, h Planck constant, in (6.64 x 10 -27 kg) helium 4 atomic mass, and c± (239 
m/s) first sound velocity in He II. Substituting with the given values above into (pQ) 
gives v c (d) = 6.42 x 10 -7 ^ m/s, as graphically shown by the broken line in Figure [TJ 
A similar formula to (pQ) was earlier derived by R. Feynman (1954) on an alternative 
theoretical argument. 

For d < 10 _5 m, v c {d) of §\§ deviates significantly from the experimental data, 
as the broken line of Figure [TJ clearly shows. Recent inelastic neutrons scattering 
experiments [31 2] have shown that, except for some layers by the wall, the bulk 
excitation spectrum of He II is not altered because of confinement even at a scale 
of ~ 70 A. So a the phonon excitation picture of He II for deriving (UJ) apparently 
retains valid for the low d region down to plausibly a few a. On the other hand, we 
observe that as d reduces, yet with d » a, the dissipation, say of energy e v( j w per 
interfacial atom, due to the fluid-wall van der Waals attraction becomes increasingly 
significant compared to the (reduced) flow energy, e s , as a result that the atomic 
population near the wall becomes increasingly a significant fraction compared to the 



1 



atomic population in the entire fluid. This is contrasted to in the wider channels 
where e v d w << £ s ; £ v dw = is one basic assumption for deriving ([1]). We below 
reexpress the critical velocity by including the VdW dissipation. 

Suppose the total VdW force acted by the wall atoms on the flow atoms as 
mapped to per interfacial-layer atom is /. The presence of the VdW force, for 
its given characteristics, should not alter the underlying QCE mechanism for the 
occurrence of a critical velocity [T]. Consider the flow has just reached its critical 
velocity v s = v c at time t = 0, and is then decelerated to v s = at t = 5t due to the 
total dissipative force on it. In this course, the dissipative work done by the VdW 
force between the wall on the flow in 5t is: 

£vdw = / fd£= m vdw v s dv s = --m vdw Vc = --(f5t)v c . ([2]) 

JO Jv c 2 2 



where f5t = m wdw v c §3§ 

is the flow momentum consumed in time 5t; m vdw is the effective fluid mass per 
interfacial-layer atom the VdW force does work on. / may be in turn expressed 
based on two phenomenological considerations: (a) / is as just said the sum of the 
VdW attractions of all of the wall atoms with all of the fluid atoms N s , mapped 
on to per collision atom out of a total Na collision atoms, thus / cx N s /Na- Or 
equivalently, / oc N s m/NA = M s /Na = m e{ , where 

77t ef = (N s /N A )m = md/Aa (gD 

(b) From general condensed matter theory we know that, / oc 1/r 6 for VdW attrac- 
tion. This function is short ranged; so atoms far away from the wall, say beyond a 
characteristic distance, do, will not contribute to /. And for d < do, the smaller the 
width d, the greater portion of the flow atoms will be subject to the VdW attraction 
from the wall. Hence, / oc (do/d) n ; u > 0. The exponent u reflects the net effect of 
the influences of (a)-(b). (do/d) u scales m cf of item (a) as 

f-(5t/v c ) = mjd /d) u . © 

The factor (5t/v c ) is multiplied to the left-hand side of © to yield a correct dimen- 
sion as according to ([3]). The equality of the two sides of ([5]) will be ensured by a 
realistic choice of the values of do and u below. Using ([5]), ([2j) writes: 

\ 2 . ©' 



Comparing (J2])' and ([2]) we have m vdw = m i 




If / = 0, then a given external energy supply would drive the flow into a velocity 
v' s = v' c at time t = 0, with a reduced critical flow energy e c J = \m ci v'^ . But now 
assuming / ^ 0, the same external energy will thus only drive the flow velocity to 
v s = v c and a reduced critical flow energy 
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Figure 1: Critical velocity v c {d) versus channel width d for superfluid He II. Circles 
represent experimental data at about 1.4 K (compiled by J. Wilks, The Properties of 
Liquid and Solid Helium, Clarendon Press, Oxford, 1967, page 391, Table 1.) Broken 
line shows QCE theoretical result (pQ) including phonon excitations only. Solid line 
shows QCE theoretical result ([8]) including a fluid-wall van der Waals attraction in 
addition to the basic phonon excitations. The solid line coincides with the broken 
line in the high d end, implying a negligible van der Waals attraction here. 
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The two flow energies are related by: 

1/2 1 2 

-m e{ v c - e vdw = -m ef v c . 
Swapping e v d w to the right side, substituting (J2])' for it, and reorganizing, we 

get: 

4' = \rn ci v' c 2 = \m c{ (l + {^pj , c 2 = Im> c 2 . ©' 

Where, 

represents a new effective mass and has a clear physical meaning. Namely, due to 
the VdW attraction from the wall, the fluid has a larger effective mass m' ; m! 

' ef ' ef 

increases with a decreasing d. 
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The threshold condition for the translational to thermal energy conversion for 
per collision atom is[H[2]: e£(= £s m in) = e ph> where e^ h = Qhc\/d. Replacing now e c s 
of the equation above by the appropriate ef, we get the threshold condition for the 
fluid translational to thermal energy conversion per collision atom in the presence 
of a fluid-wall VdW attraction: 

e' s c = e c ph . © 

Substituting ©' for e' s c and Qhc\/d for £p h as just given into ([7]), rearranging, 
we obtain the corresponding critical velocity: 




do and u are adjustable parameters and may be fixed by requiring v c (d) to agree 
with the experimental critical velocity data. The results from a least-squares fit 
(solid line, Figure H]) are do = 1.11 x 10~ 5 m and u = 1.48. Observe that, the fitted 
values of do and u affect only the bending extent of the curve v c (d) vs. d in the low 
d end, which is the characteristic region of item (ii) below. 



v c (d) divides according to the behavior of 1/y 1 + (jf^ of ([8]) in three charac- 
teristic regions: 



(i) . d » d , thus 1/y 1 + — 1- Here © identifies with CD; particularly 

notice that here v c (d) of (jHJ) is not affected by the values of do and u. The above in 
Figure [1] corresponds to in the high d end the solid and broken lines coincide. 

(ii) . d « d , thus 1/W1 + (fyY ~ (d/d ) u/2 « 1. © reduces to: 



'48/icia 1 2.98 x 1(T 3 . _, 

m/s. Q8a.) 



For the parameterization, values for do an d u are as just given from the fitting, and 
for the other quantities are as given earlier for parameterizing ([1]). The v c (d) values 
of ©,) are lowered than given by (H|) for d < do = 1.1 x 10~ 5 m, owing to the larger 
effective fluid mass here produced by the fluid-wall VdW attraction. This in Figure 
Q] corresponds to the solid line bending into a flatter one for d < 10~ 5 m. 

(iii). d ~ do- Here v c (d) undergoes a gradual switch from the steeper function 
of (HJ) to the flatter one of ([Hh.). As we see in Figure [lj do = 1.1 x 10~ 5 m represents 
just the channel width where the steeper broken line starts to give way to the flatter 
sector of the solid line. 
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